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Abstract. - We investigate geometric phase of fermion states under relative vibrations of two 
sublattices in graphene by solving time-dependent Schodinger equation using Floquet scheme. In 
a period of vibration the fermions acquire different geometric phases depending on their momenta. 
There are two regions in the momentum space; the adiabatic region where the geometric phase can 
be approximated by the Berry phase and the chaotic region where the geometric phase drastically 
fluctuates in changing parameters. The energy of fermions due to vibrations shows spikes in the 
chaotic region. The results suggest a possible dephasing mechanism which may cause classical-like 
transport properties in graphene. 



Due to advances of material science, the graphene, as a two-dimensional (2D) system 
with one layer of carbon atoms, has been fabricated recently. Graphene exhibits striking 
properties which attracted much attention of both experimentalists and theorists [1-14]. 
Particularly, the Dirac dispersion relation of electrons in graphene and the Fermi level near 
the Dirac point lead to specific features different from those in usual metals and semiconduc- 
tors. There exists a universal maximal resistivity, independent of their shapes and mobility. 
Moreover, it is found that the weak localization is strongly suppressed [6] that could be at- 
tributed to a dephasing effect similar to the phase uncertainty caused by random magnetic 
field. There are several theoretical studies using different methods addressing the unusual 
transport properties in graphene [7-14]. 

The effect of the Berry phase acquired by fermions in cyclic motions around the Dirac 
point has been noticed both theoretically [15] and experimentally [3]. The non-adiabatic 
effects and Kohn anomaly due to electron-phonon interaction have also been addressed 
recently [16]. In this paper we investigate the geometric phase of fermions induced by 
periodic vibrations between two sublattices. It is shown that in a period of vibration the 
geometric phase acquired by fermions can be expressed by the adiabatic theory only in the 
momentum region where the level splitting of fermions in the evolution path is much larger 
than huj with lo being the frequency. Outside this region the adiabatic condition is not 
satisfied. We show that the geometric phase, the fermion average energy, and the energy 
uncertainty in the non-adiabatic region are much different from those in the adiabatic region. 
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The fermion band in graphene can be described by a tight-binding Hamihonian with one 
TT orbital per site on a 2D honeycomb lattice [17]: 

H = tnn'{alfin' + a^'fln), (1) 

(nn') 

where (a„) creates (annihilates) an electron at site n, (. . .) denotes the nearest-neighbor 
(NN) sites, and t„„' is the NN hopping. The spin indices are not explicitly included. On a 
honeycomb lattice there are two sublattices, labeled as A and B. On a perfect lattice, all the 
NN hopping integrals are the same, t„„' = t^. When there is a displacement d between two 
sublattices, the lengths of three NN bonds connected to a site of sublatticc A become l'^^^ = 
yj[dsuv{e - a(-''))]2 + [Iq - dcos{0 - a^-''))]^, where j = 1, 2, 3 labels the three NN bonds, 
is the angle of d related to the x axis, lo is the original bond length, and = 0, — ^ 
for j = 1,2,3, respectively, is the original azimuth angle of the jth bond. For a small 
displacement, we can keep only the first-order terms of d and obtain /^^' ~ lo — dcos{6~a'^^'>). 
The corresponding hopping integrals become i*^-'^(d) = to + \dcos{6 — a^^'), with A being 
the coefficient of linear dependence of i^^^(d) on /^-'X 

By using Bloch transformation for electron operators the original Hamiltonian becomes 

Hq ~ ^ to 

k 

where a(6)k is annihilation operator of electron on sublatticc A(B) with Bloch wavcvector 
k. Expanding it to the first order of k around two irreducible K points ^0,^ ^^^ ^ in the 
Brillouin zone, one has 

k 

where 1 and fz are unit and Pauli matrices acting on two valleys at two irreducible K 
points, and a^^y are Pauli matrices on two sublattices. The interaction between electrons 
and lattice displacement d is 

-f^i = ^ (t„„'(d) - to)(aJ,a„' a|j,a„), (4) 

where t„„' (d) is the hopping integral between sites n and r\! under displacement d. Keeping 
only the terms of the first order of d, the interaction Hamiltonian becomes 

Hi = —'^^(dji® a^c^ dyT^® Oy). (5) 

k 

It is similar to the expression previously derived based on a valence- force-field model [18]. 
We only consider in-plane optical modes at the long wavelength limit, as the modes with 
short wave lengths induce large momentum transfer which causes the electron states out 
of the Dirac fermion region, the vertical optical modes only quadratically couple with the 
fermions, and the coupling between acoustic modes and electrons is even more weak as the 
nearest-neighbor bond lengths are almost unchanged. We denote the in-plane vibrations of 
the relative coordinates as = ^ cos(a;t -I- ax) and dy = ^ cos(tjt -I- ay), where ^^^^ and 
ax(y) are amplitude and initial phase of the vibration in the x{y) direction, respectively. 
The Hamiltonian for electrons becomes time-dependent: 

H{t) = 2 5Z [i'^olokyfz + rjx cos(wt -f ax)i) <E) (Jx - {hlokxi - % cos(wt + ay)Tz) <E) (Jy] ■ 

k 

(6) 
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■y'o 



at^k 



H.c. 
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Fig. 1: (Color online) E\ ^ as functions of d for different values of k. Green circle shows the 
evolution path in the vibrations. 

For a given valley labeled by + or — , the Hamiltonian is: 

H±{t) = ^ X] {^o^olifcyCT^ - kxUy) + ^x. cos{ujt + ax)cTx ± cosiiut + ay)ay} . (7) 

k 

At the adiabatic limit, the instantaneous eigenenergies at time t can be obtained by diag- 
onalizing H±{t): 

E^^\i) = i-lT' -\tM±ky + '^k^)+VxCOs{ujt + ax)T'^VvCOs{ujt + ay)\, (8) 

where m = \ and m = 2 correspond to bands below and above the Dirac point, respectively. 
For given (fc^:, ky)^ these instantaneous eigenenergies exhibit conic dependence on parameters 
dx and dy, as shown in Fig. 1. The diabolical point, where the poles of two cones coincide, 
is determined by {dx = ^tokylo/ X, dy = zLtokxlo / X) , where plus and minus signs refer to 
two valleys. The vibrations of dx and dy result in circular motion with an elliptic track 
in the dx ~ dy plane. From the Berry theorem, a cyclic motion along a close track in a 
2D parameter space can cause a Berry phase of ±7r in a wavcfunction whose instantaneous 
eigenenergy has a diabolical point enclosed in this circle [19]. Since the position of the 
diabolic point in the dx — dy plane is determined by {kx,ky), the Berry phase acquired by 
the fermions with different momenta are different: It is ±7r for the fcrmions whose diabolic 
point {dx = ^takylo/ X, dy = ±tQkxlo/X) is enclosed in the track, while it is zero for the 
states with the diabolic point in the dx — dy parameter space outside the track. 

The above conclusion about the Berry phase is valid only in the adiabatic condition, i.e., 
the energy difference between two bands is always much larger than fiuj in the track. There 
exist fermion states for which the adiabatic condition can not be satisfied. In this case we 
have to investigate the geometric phase from solutions of time-dependent Hamiltonian. In 
a periodic time-dependent Hamiltonians the wavefunctions can be expressed in the Floquct 
form [20,21]. This scheme has been used to solve the states in a laser field in calculation of 
the harmonics generation in graphene [22] and in nanotube [23]. In periodic vibrations of 
dx and dy the electron wavefunction can be written as 

^^{l,k,t) = e-''^^*u^{l,k,t), (9) 
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where I is an index of electron states for given k, and u^{l, k, t) is a periodic function of t 
with period T = Performing Fourier transformation for {I, k, t) with respect to t, one 
has 

.±n u f ^ ^_i„„t 



--il,k,t) 



(10) 



where and i?^ are components on two sublatticcs. Substituting the Floquet wave 
function into the time-dependent Schrodingcr equation [H±{t) — ifidt]4'^ (lykjt) = 0, we 
obtain a set of Unear homogeneous equations for the components 



2 



{±ky+ik^)Bl 



4^ 



(±fc,-ifc,)4± + -(e-==,7,±ie-«,7y)^t 



, + -(e--^r,,±ie"-«,7y)4^ 



+1 



(e±+7i?iw)4±, 
(11) 

(e± + nritj)S±. 
. 

From the requirement of nonzero solutions, for given k one can solve the discrete quasicner- 
gies e*(/,k) and the corresponding Floquet states ■;/'*(/, k, t). It is noteworthy that two 
Floquet states whose quasienergies differ by nhuj with n being an integer are physically 
equivalent states [21]. So the quasienergies of physically different Floquet states can be 
reduced into region - ^ <£=*=< 

In a period the wavefunction ?/;*(/, k,t ~ T) ~ e*"^ ^''''V*(^k, i = 0) acquires a phase 
0±(Z,k). It consists of two parts: 0±(?,k) = a± (/, k) +/3±(/, k), where a±(/,k) = ^E^^l^ 

is the dynamical phase with i?±(Z,k) = ^ /J" di('0=^(/, k, t)|if± (i)|V'=^ (/, k, t)), and /3± (/, k) = 
(/)*(Z,k) — a*(Z,k) is the geometrical phase [24]. Under the periodic vibrations the Flo- 
quet state V'*(^k, t) is stationary one which returns to the initial state other than phase 

(/)^(l, k) = ~^^^^in7^ after a period of evolution. At the same time, the energy is no longer 
a good quantum number and the average energy of stationary state '0^(/,k, t) can be cal- 
culated as 



^=^(/,k) = e^{l,'k) + h}_^nw{\A^{l,\i)\^ 

n 

From this one obtains the geometric phase 



/3±(Z,k) = 27r^n(|^±(Z,k)p + |i?±(Z,k) 



(13) 



(14) 



The energy uncertainty caused by the vibration can be specified by the standard variance 

A±(Z,k)s(i?2)-(i?)2 



^n(|Aj(;,k)p + |B±(;,k)|2) 



(15) 



From the superposition principle any time-dependent state can be expressed as a linear 
combination of the Floquet states. But this linear combination is usually not a stationary 
state, i.e., the state could not return to its initial one with only a phase difference after 
a period. For such states the geometric phase can not be defined. So in this paper we 
only consider the geometric phase for the Floquet states. Since for given valley and given 
momentum there are only two unknowns having the same phonon number n in Eqs. (11) 
and (12), the number of quasienergies within range of [— is 2, corresponding to the 
lower and upper bands of the Dirac fermions with labels I = 1 and 1 = 2, respectively. 

Now we begin to investigate the properties of stationary states in the periodic vibrations. 
We are interested in: (i) the geometric phases acquired by fermion states with various 
momenta in a period of vibrations; (ii) the deviation of the average energy from the Dirac 
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Fig. 2: (Color online) Geometric phase acquired in a period of vibrations by Dirac fermions in 
lower and upper bands of valley "+" as a function of the fermion momentum. The parameters for 
vibrations are; rjx ~ riy = rjo, Ox = 0, and Oy — kx = and units of ky are set to be -j^fl- so 
that the border between the ±7r Berry phase and the zero Berry phase is at fcy = 1. The inset of 
(c) shows the width of the chaotic region Sky as a function of the frequency. 



dispersion relation and the energy uncertainty caused by the vibration. From (i) we can 
shed some light on the phase variation of Dirac fermions. From (ii) we can see the essential 
effect of vibrations on the basic energy spectrum. 

In Fig. 2 we plot the geometric phase acquired by fermion states in a period of vibrations 
versus the fermion momentum. There is a ^ phase difference between the vibrations in 
the X and y directions and their amplitudes are the same, i.e., r/x = Vy = Vo- So this 
is a circular vibration of the relative coordinates of two sublattices as shown in Fig. 1. 
According to the Berry theorem, the adiabatic Berry phase is ±7r if |k| < and is zero 
if |k| > We note that from the calculations beyond the adiabatic approximation the 

obtained geometric phase is still roughly equal to the Berry phase: for |k| <^ it is nearly 
±7r and for |k| 3> -pj- it is almost zero. However, now there appears a "chaotic" region 
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around the border |k| = where the geometric phase randomly oscihates in changing 
the momentum. This is a direct consequence of the non-adiabaticity, as in this region the 
adiabatic condition, that the energy difference between two bands is much larger than fiio 
in the evolution path, is not satisfied. As can be seen from the inset of Fig. 2(c), the 
width of this chaotic region is increased by increasing the frequency. The chaotic nature 
of the geometric phase reflects the phase uncertainty of the single-fermion states in this 
momentum region during vibrations. Such vibrations may exist in the graphene due to 
thermal excitations, or due to the zero-point fluctuations. Especially, the Berry phase effect 
may produce Born-Huang vibronic centrifugal term (see below) which may cause nonzero 
vibrations even at very low temperature. So such a fermion-lattice interaction can play 
roles of a dephasing mechanism for the single-fermion states. We also note that in the whole 
region, even including the chaotic region, the geometric phases acquired by fermions of the 
same momentum in two bands are opposite to each other. This phase compensation effect 
of two bands implies that this dephasing mechanism has no effect if the fermions in two 
bands are paired. The exchange of two bands leads to the sign reverse of the geometric 
phase. From the symmetry shown in Fig. f we can see that an equivalent operator which 
can also cause the sign reverse of the geometric phase is the reverse of the direction of the 
circular vibration from clockwise to counterclockwise or vice versa. 

Another significant effect of vibrations is the change of fermion energies. As the single- 
particle energy is not a good quantum number, in Fig. 3 we plot the average energy of 
single-fermion states as a function of the momentum. From the comparison with the Dirac 
dispersion relation and with Fig. 2, the energy difference between two bands is enlarged near 
the Dirac point and for |k| ^ where the geometric phase is near its adiabatic value, 
±7r or zero, while this energy difference shrinks in the chaotic region where the adiabatic 
approximation can not be used. Such opposite behaviors in these two regions imply that the 
failure of the adiabatic theory in the chaotic region has much more profound meaning than 
reflected from the values of geometric phase. Physically, in the chaotic region the oscillations 
of electrons between two bands are not able to follow the vibrations, this causes the loss of 
distinguishability of the two bands and the resultant states trend to take average energies 
in between, leading to a smaller energy spacing. On the contrary, in a quantum description 
the coupling of two states always enlarges their energy spacing. So one may expect that the 
chaotic region corresponds to a classical-like behavior and has a maximum uncertainty of 
single-fermion energies. 

The total average energy of fermions Ve is then also changed and becomes dependent on 
the vibration amplitude 770- In the inset of Fig. 3(a) we plot the effective potential V{r]o) = 
^9Vo + ^eivo) where the first term is the ordinary elastic potential. The minimum is shifted 
from the center (Born- Huang vibronic centrifugal term) , leading to nonzero vibrations even 
at very low temperature. 

We calculate the standard variance of Eq. (15) and show the results in Fig. 4. As 
expected, the standard variance exhibits spikes in the chaotic region, and the energy un- 
certainty increases by increasing the amplitudes of vibrations. Except for the spikes in the 
chaotic region, the standard variance globally increases by increasing the momentum. Es- 
pecially, the energy uncertainty becomes zero at the Dirac point for all the investigated 
vibration amplitudes. This originates from the fact that the the energy uncertainty is due 
to the oscillation of the fermion states between two bands, but at the Dirac point the two 
bands coincide, leading to a zero amplitude of the oscillation. This structure, including the 
spikes and the globally increasing background with the momentum, may catch major fea- 
tures of linewidth distribution obtained from the angle-resolved photoemission spectroscopy 
in graphene [25]. Spikes also appear in the average energy of electrons in quasi-ID nanotube 
under periodic laser field [23] , but in the present case the spikes occur in 2D momentum space 
originating from specific electron-phonon interaction in graphene and are closely related to 
the geometric phase. 

The evolution path plays a crucial role on the distribution of adiabatic and chaotic 
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Fig. 3: (Color online) Average energy of fermions in lower and upper bands of valley "+" as a 
function of the fermion momentum. The parameters for vibrations are rjx = rjy = rjo , ax = 0, and 
ay = kx = 0. The blue dotted lines represent eigenenergies without vibrations. Inset in (a): 
The effective potential of vibrations as a function of rescaled displacement 770 for different carrier 
densities. For curves from bottom to top the extra electron densities above the Dirac point are: 0, 
0.00054, 0.00108, 0.00162, 0.00216, 0.0027 per atom. Hlj = 0.0037to and g = 0.25/fo- 



regions in the momentum space. To show this in Fig. 5 we plot the geometric phase as a 
function of and ky for an eUiptic path. The chaotic region forms an orbicular area with 

a nearly fixed width along a loop defined by kx — ^.nd ky ~ ^^^T^I^- The width 

depends on uj as shown in the inset of Fig. 2(c). As a result, the adiabatic region with a ±7r 
Berry phase is compressed by reducing the short axis of the ellipse, but the chaotic region 
with random geometric phases can exist even for a linear vibration. 

In summary, we investigate the geometric phase and the change of dispersion relation 
of fermion states under relative vibrations of two sublattices in graphene. For a circular 
vibration there are two regions in the momentum space: in one region the geometric phase 
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Fig. 4: (Color online) Standard variance of single-fermion energy of valley "+" as a function of 
the fermion momentum. The parameters for vibrations are rj^ ~ rjy — rjo, ~ 0, and Oy = 
kx = 0. The standard variance is the same for lower and upper bands. 



acquired by the fermion states can be approximated by the adiabatic Berry phase and the 
average energy spacing between two bands is enlarged; in the other region the geometric 
phase exhibits random oscillations in changing the momentum, much different from the 
adiabatic value, and the average energy spacing shrinks. The energy uncertainty of fermions 
due to vibrations shows spikes in the chaotic region. For elliptic paths of vibrations the 
distribution of the adiabatic and chaotic regions crucially depends on the path shapes. 
The results suggest a possible dephasing mechanism which causes classical-like transport 
properties in graphene. 
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